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Abstract 

We report a measurement of 4> meson production in central Pb+Au collisions at 
Ejab/A=158 GeV. For the first time in heavy-ion collisions, </> mesons were reconstructed in the same 
experiment both in the K + K _ and the dilepton decay channel. Near mid-rapidity, this yields ra- 
pidity densities, corrected for production at the same rapidity value, of 2.05±0.14(stat)±0.25(syst) 
and 2.04±0.49(stat)±0.32(syst), respectively. The shape of the measured transverse momentum 
spectra is also in close agreement in both decay channels. The data rule out a possible enhancement 
of the 4> yield in the leptonic over the hadronic channel by a factor larger than 1.6 at 95% CL. 

PACS numbers: 25.75.-q, 25.75.Dw, 13.85.Ni, 13.85.Qk 
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In ultrarelativistic nucleus-nucleus collisions matter is created at high temperature and 
density. For top SPS energies and higher, the corresponding fireball is very likely Q in a 
deconfined state as predicted by solutions of Quantum Chromodynamics (QCD) on a space- 
time lattice P| • Detailed studies of hadron abundances imply [3j that their yield is frozen very 
close to the predicted phase boundary between confined and deconfined matter. Strangeness 
enhancement has been suggested 4J as a signature of a deconfined stage. In this context 
it is important to understand the production of mesons which carry hidden strangeness. 
Furthermore, near the phase boundary the d> meson mass, width and its branching ratios 
into kaons and leptons might be modified 0] . Due to final state interactions of the kaons, 
such a modification would probably only be visible in the lepton decay channel, which 
could result in different <fi yields depending on which decay channel is studied |7|. Indeed, 
the experimental situation concerning <fr production at the top Super Proton Synchrotron 
(SPS) energy provides hints for such a scenario as the NA50 collaboration reported a <fi 
yield measured via dileptons [s| which exceeded the yield determined by NA49 in the K + K~ 
channel Q] by factors between two and four ^] in the common transverse momentum range. 
Further, the m t spectra exhibit a different inverse slope parameter [9J], 305 ±15 MeV in 
NA49 and 218 ± 6 MeV in NA50, fitted in their m< acceptance regions. 



The upgrade of the CERES experiment 
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14j makes possible for the first time 



to simultaneously study the leptonic and the charged kaon decay modes of the (f) meson at 
the SPS, thus shedding light onto the (f> puzzle. In this paper, we present results of 4> mesons 
reconstructed both in the charged kaon (K + K~) and in the dilepton (e + e~) decay mode. 
The kaon (dilepton) analysis uses 24 (18) million 158 AGeV/c Pb on Au collisions taken at 
the most central 7% of the geometrical cross section. 

The CERES experiment is optimized to measure low mass electron pairs close to mid- 



rapidity (2.1< 77 <2.65) with full azimuthal coverage A vertex tele- 

scope, composed of two Silicon Drift Detectors (SDD) positioned at 10 cm and 13.8 cm 
downstream of a segmented Au target (thickness equivalent to 1.33% interaction length), 
provides a precise vertex reconstruction. Two Ring Imaging CHerenkov (RICH) detectors, 
operated at a high threshold (7^=32), provide electron identification in a large hadronic 
background. The new radial-drift TPC, positioned downstream of the original spectrome- 
ter inside a magnetic field provides the momentum of charged particles with a resolution 
Ap/p ~ ((2%) 2 + (1% • ^(GeV/c)) 2 ) 1 / 2 and additional electron identification via dE/dx. 
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Charged particles from the target are reconstructed by matching track segments in the SDD 
and in the TPC using a momentum-dependent matching window. Matching to a ring in the 
RICH is also required for the dilepton analysis. 

To study the meson in the charged kaon (K + K~) decay mode, all charged particles 
get assigned the kaon mass (no particle identification is used). Only a conservative upper 
cut in the dE/dx signal (corresponding to 90% of the Fermi plateau value) for momenta 
between 1.25 GeV/c and 4 GeV/c, suppressing 83% of the electrons, is applied to enhance 
the kaon content of the sample. Tracks in the TPC are required to contain more than 12 hits 
of maximum 20 possible to ensure good momentum resolution. Tracks in the geometrical 
acceptance 0.13 rad < 9 < 0.24 rad with a transverse momentum p t larger than 0.25 GeV/c 
are selected. To reduce the contamination from other particle species, cuts in the Podolanski- 
Armenteros parameter j^cj and in the opening angle between the kaons are applied 21 1. 



The meson in the kaon decay mode is studied in the rapidity interval 2.0 < y^ < 2.4 for 
pf > 0.75 GeV/c. 

To obtain the <fi transverse momentum spectrum the invariant mass distributions of K + K~ 
pairs were accumulated in (pf,y^) bins. The invariant mass distributions of the combinatorial 
background are calculated using the mixed-event technique for each (pt,y^) bin. One such 
invariant mass distribution, after background subtraction, is presented in Fig. ^ Residual 
background appears in the low pf bins. The yield of the mesons is, therefore, determined 
by fitting a relativistic Breit-Wigner distribution with parameters taken from the Particle 
Data Group compilation [3] (convoluted with the experimental resolution function obtained 
by a Monte Carlo simulation) superimposed on a linear background, to the measured line 
shape. The signal to background ratios vary from 1/2000 to 1/180 with increasing pf. The 
signal is integrated in the mass range between 1.0 GeV/c 2 and 1.05 GeV/c 2 . The resulting 
yields need to be corrected for acceptance and efficiency. 

To study the meson in the dilepton (e + e~) decay mode, electrons are identified among 
all charged hadrons using the RICH detectors and the dE/dx signal in the TPC. The two 
RICH detectors are used in a combined mode. Cherenkov rings with asymptotic radius are 
identified using a Hough transformation (see sect. 3.2.4 in ref. ^j]). A pion rejection factor 



of 2000 is achieved in the RICH |23j for an electron efficiency of 0.70 with the quality cuts 
applied in the analysis (24, [2^. The TPC electron selection is done based on the dE/dx 
signal and its resolution. The combined pion rejection factor varies from 4xl0 4 to 1.8 xlO 4 
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FIG. 1: Invariant-mass spectra of K + K pairs after background subtraction for the (pf,y^) bin 
1.5 GeV/c < pf < 1.75 GeV/c and 2.2 < < 2.4. 

or momenta between 1 GeV/c and 2.5 GeV/c for total electron efficiencies of 68% and 66% 



24|, respectively. 



The main difficulties of reconstructing mesons in the di-electron channel are the low 
branching ratio of this decay mode and the large amount of combinatorial background from 
7 conversions and Dalitz decays of neutral mesons. Therefore, besides a very good electron 
identification, removal of electron pairs from 7 conversions and tt° Dalitz decays from the 
sample is vital to reduce the combinatorial background. The details of the rejection strategy 



are explained in 
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24, 



25 



271. Electron tracks in the geometrical acceptance 0.14 rad < 

nn 



9 < 0.243 rad and with a transverse momentum p t > 0.2 GeV/c are selected 24j,|25j. The 
meson in the dilepton decay mode is studied in the rapidity interval 2.1 < < 2.65. The 
signal in the e + e~ channel is obtained from the invariant-mass distributions of unlike-sign 
pairs after full rejection and subtraction of the mixed-event combinatorial background. The 
mixed-event background is normalized to the like-sign pair background in the mass region 
0.2 GeV/c 2 < m e + e - <1.6 GeV/c 2 . The integrated yield in the mass range between 0.9 and 
1.1 GeV/c 2 is 229±53 with a signal to background ratio of 1/12 and needs to be corrected 
for acceptance, reconstruction efficiency and physics background under the peak. 

In order to determine the acceptance, decay and efficiency corrections, Geant [26] sim- 
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ulations containing the description of the CERES experiment are used. The Monte Carlo 
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Sources of 6 mesons with 



simulation is tuned to reproduce all aspects of the data 
realistic transverse momentum and rapidity distributions were embedded into real events to 
simulate the background. The simulation of the dilepton channel is done using the GENE- 
SIS simulation code |28j]. All (pt,y^) bins with an acceptance larger than 0.5 are used in the 
K + K~ analysis. The reconstruction efficiency is of the order of 0.4 with slight dependence 
on the transverse momentum, given by the decay in flight of the charged kaons, the pair cuts 
applied and the single track efficiency. The acceptance for the e + e~ pairs varies from 0.165 
to 0.185 in the measured rapidity range and the reconstruction efficiency after full rejection 
varies from 0.22 to 0.145 with the transverse momentum. An additional loss of 25% is due 
to bremsstrahlung events that lead to a reduction of the mass by more than 100 MeV. 

The e + e~ invariant-mass spectrum corrected for efficiency and normalized to the number 
of charged particles in the acceptance (^tJ is shown in Fig. E] together with the expectations 
from the hadron decay cocktail [28]. The meson yield in the e + e~ channel is determined by 
integrating the invariant mass spectra in the mass region between 0.9 and 1.1 GeV/ c 2 in three 
transverse momentum bins. The p meson could extend into this mass range if its spectral 
function is modified in the medium. Dileptons from the QGP phase also contribute to the 
physics background in this mass range. The sum of these two contributions is estimated to 
be 35% of the total yield in this mass region by inspecting theoretical models that include 
in-medium s pre ading of the p width due to 27r processes and the dilepton yield from the 
QGP phase |3j|. If processes involving 4 or 6 pions would be included in |3Jj the physics 
background could be larger. Using another model j^] that also describes the CERES data 
gives a physics background contribution of 37%. The measured yield has been scaled by 
the smaller factor to correct for the physics background. The charm contribution that is 
smaller than 3% has been neglected. 

The efficiency- and acceptance-corrected meson yield is shown in Fig. El for both decay 
modes as a function of transverse momentum. The error bars shown are statistical. The 
systematic errors in the charged kaon analysis (plotted as boxes) contain contributions from 
variation of the pair cuts, from the difference between the integral from the fit and from data 
points in the given interval, from the dependence of the yield on the fit window, and the 
function used to fit the residual background. The systematic errors in the dilepton analysis 
are ±16%. They contain the contribution from the analysis cuts, background substraction 
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FIG. 2: Invariant mass spectrum of e + e~ pairs per charged particle compared to the hadron decay 
cocktail (thick solid line). Also plotted is a model calculation assuming in medium spread p width 
plus the dilepton yield from the QGP phase (dashed line) used to extract the physics background 
in the <j) mass region. 

and variation of mass range. They are not drawn because the errors in the dilepton channel 
are totally dominated by statistics. 

When the spectra are fitted with the function 

d N/dMy = T-(T + m,) ' Pt ' 6XP ( — ] 

an inverse slope parameter of T = 273 ± 9(stat)±10(sys) MeV and a rapidity density dN/dy 
of 2.05 ± 0.14(stat)± 0.25(sys) in the K + K~ decay mode and T = 306 ± 82 MeV and 
dN/dy = 2.19 ± 0.52(stat)±0.34(syst) in the dilepton decay mode are obtained averaging 
over their corresponding rapidity intervals. In order to compare the rapidity densities in the 
two decay channels the dilepton data are scaled by 10] (dA^/dy) 2-2 - 4 =0.93-(d^/dy) 2 - 1 ~ 2 - 65 . 
The meson yields and inverse slope parameters obtained in both decay modes agree within 
the errors. A <fi yield in the e + e~ decay mode larger than 1.6 times the yield on the K + K 
decay mode is excluded at 95% CL (statistical and systematic errors in both decay channels 
added in quadrature). 

Q 

Moreover, the CERES results can be compared to the existing Pb-Pb systematics [9] after 
accounting for the different measurement conditions. The NA49 measurement was done at 
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FIG. 3: Transverse momentum spectrum of (j) mesons corrected for acceptance and efficiency 
reconstructed in the e + e~ decay mode (closed symbols) and in the K + K~ decay channels (open 
symbols). The spectrum in the e + e~ decay mode has been scaled by 0.93 before comparing to the 
K + K~ decay channels to account for the different rapidity intervals covered in the two analyses. 
Systematic errors in the kaon decay channel are shown as boxes. 



4% centrality and covered a rapidity range from 3 to 3.8 units [10]. We use the measurec 
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increase of charged hadron (h~) multiplicity with centrality ( h^ % /h^ % = 1.092±0.014 
to correct the <fi meson yield for the different centralities between the two experiments. From 
the NA49 <b meson rapidity distribution a ratio of (dN /dy) 3 - 3 - 8 /(dN /dy) 2 -- 2 - 4 = 1.07 ± 
0.11 is obtained to account for the different rapidity coverage. Thus, a global scaling factor 
of 1.17 ± 0.12 is applied to the combined CERES data of Fig. [31 to make the comparison to 
the systematics of ^ . In Fig. 0] the transverse mass spectrum obtained by CERES after 
scaling is plotted together with the NA49 and NA50 data. The meson yields agree within 
the errors with the NA49 results. So does the yield in the K + K~ extrapolated down to 
Pt=0 using the measured inverse slope parameter. On the other hand, CERES data in the 
K + K~ channel do not agree with NA50 results in the common p t region. This experiment 
however measures the leptonic channel. The extrapolation of NA50 results down to the 
region where CERES measures the dilepton channel does not agree either. As stated above, 
in the CERES measurement the two decay modes agree. Possible differences of maximum 
40-50% as expected by models including only rescattering of the kaons 7} or of maximum 
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FIG. 4: Transverse mass distribution of (j) mesons measured in the charged kaon (triangles) and 
in the dilepton (circles) decay mode after scaling (see text) compared to the results from NA49 
(squares) and NA50 (diamons). Error bars are statistical and systematic errors in the kaon channel 
are shown as boxes. 



70% at the lowest p t (pt < 0.3 GeV/c) expected by models including medium modifications 
of the <fi mesons and kaons like the AMPT model j^J cannot be ruled out by the CERES 
results. 

To conclude, meson production has been measured simultaneously in both decay chan- 
nels for the first time at the CERN SPS. The yield and the inverse slope parameter in both 
decay modes agree within the errors. Our results are in agreement with the results from 
NA49 measured in the kaon channel. A yield in the e + e~ decay mode larger than 1.6 times 
the yield on the K + K~ is excluded at 95% CL, therefore the large discrepancy observed 
previously is not observed in the CERES data. The theoretical predictions in ?} 3^ are 
consistent with our data but with our precision we cannot distinguish between the different 
effects. 

The CERES collaboration acknowledges the good performance of the CERN PS and 
SPS accelerators as well as the support from the EST division. We would like to thank 
R. Campagnolo, L. Musa, A. Przybyla, W. Seipp and B. Windelband for their contribution 
during construction and commissioning of the TPC and during data taking. We are grateful 



9 



for excellent support by the CERN IT division for the central data recording and data 
proccesing. This work was supported by GSI, Darmstadt, the German BMBF, the German 
VH-VI 146, the US DoE, the Israeli Science Foundation, and the MINERVA Foundation. 



[1] See, e.g., Proc. Quark Matter 2004, J. Phys. G: Nucl. Part. Phys. 30. 

[2] F. Karsch, J. Phys. G: Nucl. Part. Phys. 31 (2005) S633. F. Karsch, Nucl. Phys. A 698 (2002) 
199c. 

[3] P. Braun-Munzinger, J. Stachel, C. Wetterich, Phys. Lett. B 596 (2004) 61. 

[4] P. Koch, B. Miiller, J. Rafelski Phys. Rep. 142 (1986) 167. 

[5] D. Lissauer and E. V. Sfmryak, Phys. Lett. B 253 (1991) 15. 

[6] F. Klingl, T. Waas, W. Weise, Phys. Lett. B 431 (1998) 254. 

[7] S.C. Johnson, B.V. Jacak, A. Drees, Eur. Phys. J. C18 (2001) 645. 

[8] N. Willis et al, NA50 Collaboration, Nucl. Phys. A 661 (1999) 543c; B. Alessandro et al, 

Phys. Lett. B 555 (2003) 147. 

[9] D. R6hrich, J. Phys. G 27 (2001) 355. 

[10] S.V. Afanasiev et al, NA49 Collaboration Phys. Lett. B 491 (2000) 59. 

[11] Addendum to proposal SPSLC/P280: CERN/SPSLC 96-35/P280 Add.l. 

[12] Technical Note on the NA45/CERES upgrade. CERN/SPSLC 96-50 (1996). 

[13] A. Man'n for the CERES Collaboration, Nucl. Phys. A 661 (1999) 673c. 

[14] A. Man'n et al, CERES Collaboration, J. Phys. G: Nucl. Part. Phys. 30 (2004) S709. 

[15] G. Agakichiev et al, CERES Collaboration, Eur. Phys. J. C 4 (1998) 231. 

[16] G. Agakichiev et al, CERES Collaboration, Phys. Rev. Lett. 75 (1995) 1272. 

[17] G. Agakichiev et al, CERES Collaboration, Phys. Lett. B 422 (1998) 405; B. Lenkeit for the 
CERES Collaboration, Nucl. Phys. A 661 (1999) 23c; J. P. Wessels for the CERES Collabora- 
tion, Nucl. Phys. A 715 (2003) 607c. 

[18] G. Agakichiev et al, CERES Collaboration, Eur. Phys. J. C 41 (2005) 475. 

[19] D. Adamova et al, CERES Collaboration, Phys. Rev. Lett. 91 (2003) 042301. 

[20] J. Podolanski and R. Armenteros, Phil. Mag. 45 (1954) 13. 

[21] CERES Collaboration, Reconstruction of 4> mesons in the K + K~ channel, in preparation. 

[22] S. Eidelmal et al, Phys. Lett. B 592 (2004) 1. 

10 



[23] O. Busch, Doctoral Thesis, in preparation, TU Darmstadt (2005). 
[24] S. Yurevich, Doctoral Thesis, in preparation, University of Heidelberg (2005). 
[25] A. Cherlin, Doctoral Thesis, in preparation, Weizmann Institute of Science (2005). 
[26] GEANT, Detector Description and Simulation Tool, CERN Program Library Long Writeup 
W5013. 

[27] D. Miskowiec for the CERES Collaboration. Proc. Quark Matter 2005 Nucl. Phys. A in print, 
nucl-ex/0511010 

[28] H. Sako for the CERES Collaboration. Technical Report 03-25 (2000). 
[29] V. Friese for the NA49 Collaboration, Nucl. Phys. A 698 (2002) 487c. 
[30] F. Sikler for the NA49 Collaboration, Nucl. Phys. A 661 (1999) 45c. 

[31] R. Rapp and J. Wambach, Adv. Nucl. Phys. 25 (2000) 1; R. Rapp, private comunication; 

E. Braaten, R.D. Pisarski and T.-C. Yuan, Phys. Rev. Lett. 64 (1990) 2242. 
[32] K. Gallmeister et al., Nucl. Phys. A 688 (2001) 939; and B. Kampfer, private communication. 
[33] P. Braun-Munzinger, D. Miskowiec, A. Drees, C. Lourengo, Eur. Phys. J. C 1 (1998) 123. 
[34] S. Pal, CM. Ko, Zi-wei Lin, Nucl. Phys. A 707 (2002) 529. 



11 



